Magnetic structure of the edge-sharing copper oxide chain compound NaCu202 
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Single-crystal neutron diffraction has been used to determine the incommensurate magnetic struc- 
ture of NaCu202, a compound built up of chains of edge-sharing Cu04 plaquettes. Magnetic struc- 
tures compatible with the lattice symmetry were identified by a group-theoretical analysis, and 
their magnetic structure factors were compared to the experimentally observed Bragg intensities. In 
conjunction with other experimental data, this analysis yields an elliptical helix structure in which 
both the helicity and the polarization plane alternate among copper-oxide chains. This magnetic 
ground state is discussed in the context of the recently reported multiferroic properties of other 
copper-oxide chain compounds. 



PACS numbers: 75.25.-j, 75.85. -l-t, 75.50.Ee, 75.10.-b 

I. INTRODUCTION 

Research on compounds with coupled spontaneous 
magnetic and electric polarization (multiferroics) has re- 
cently focused renewed attention on helical magnetism. 
A leading theory of multiferroicity predicts that helical 
magnetic order in insulators induces an electric polariza- 
tion of the form 

P cx ^ny- X (S, X Sj ), (1) 

(u) 

where fiij is a vector connecting spins Sij on exchange 
bond (ij)^^ Clearly, detailed information about the 
magnetic structure and crystal symmetry is required to 
evaluate the consequences of Eq. 1 for any given com- 
pound. The theory explains the observation of ferroelec- 
tricity in compounds such as TbMnOa and CuO, which 
exhibit helicoidal states with nonvanishing projection 
of the propagation vector onto the polarization plane.* 
While these compounds contain three-dimensional (3D) 
networks of exchange bonds, magnetic insulators with 
quasi- ID electronic structure have the potential to serve 
as particularly instructive model systems for multifer- 
roicity. The discovery of a macroscopic electric polariza- 
tion in the compounds LiCu202 (Ref. |S| and LiCuV04 
(Ref. E]), which are built up of chains of edge-sharing 
Cu04 plaquettes, has therefore generated significant at- 
tention. Because of the competition between the anoma- 
lously small nearest-neighbor exchange coupling and the 
stronger next-nearest neighbor coupling of spin-1/2 cop- 
per ions along the chains, both compounds exhibit in- 
commensurate magnetic correlations, and interchain in- 
teractions induce helicoidal 3D long-range order at low 
temperatures. In addition to multiferroicity, these mate- 
rials are also of interest in the context of research on the 
interplay between spin and charge correlations along the 
copper-oxide chains .^^ll^ 



In spite of their simple electronic structure, the mul- 
tiferroic properties of LiCu2 02 and LiCuV04 are still 
poorly understood. Both compounds are orthorhom- 
bic, with copper-oxide chains running along the 6-axis. 
Based on neutron diffraction experiments, the polariza- 
tion plane of the helix in LiCu2 02 was originally reported 
to lie within the ab plane ,'^^1 while P was found to be along 
c, in disagreement with Eq. 1. El Later polarized-neutron 
experiments^^ indicated a &c-polarized helix, which is 
consistent with the electrical polarization according to 
Eq. 1. However, the behavior of P in an applied mag- 
netic field H was found to be difficult to reconcile with 
this scenario, although interchain magneto-electric 
coupling may offer a possible solution.J^ Recent single- 
crystal neutron diffraction datcP^ have been interpreted 
as evidence of an elliptical helix state with polarization 
plane tilted by an angle of 45° with respect to the foe- 
plane. The consequences of this structure for the electric 
polarization remain to be assessed. In LiCuV04, the he- 
lix wa s repo rted to be polarized within the a6-plane for 
H = 0,' ^^l^^ l which is consistent with the observation that 
P II a. However, deviations from the theoretical predic- 
tions were again noted for H ^ O.^^^The apparent failure 
of Eq. 1 as a description of the magnetic field dependence 
of the dielectric properties of LiCu202 and LiCuV04 
has motivated an alternative scenario according to which 
the observed electric polarization is due to defects gen- 
erated by Li-Cu intersubstitution.'i^E^ While such de- 
fects are indeed q uite c ommon due to the similar ionic 
radii of Li and Cu,E2lil this scenario has been contested 
based on work on ostensibly stoichiometric LiCu202 sin- 
gle crystals.'^ The situation thus remains unresolved. 

Here we present an investigation of the magnetic struc- 
ture of NaCu202, a compound that is isostructural 
and isolelectronic to LiCu202 and, like LiCu202 and 
LiCuV04, exhibits incommensurate magnetic order at 
low temperatures-f^SHUl Unlike LiCu202 and LiCuV04, 
however, NaCu202 is not intrinsically prone to disor- 
der because of the different sizes of Na and Cu ions. 
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and Na-Cu intersubstitution can be ruled out with high 
sensitivity!^ A recent investigatiorP^l came to the conclu- 
sion that ferroelectricity is absent in NaCu202. Detailed 
information about the magnetic structure is required in 
order to assess the origin of the qualitatively different 
electric properties of LiCu2 02 and LiCuV04 on the one 
hand, and NaCu202 on the other hand. However, while 
a comprehensive set of single-crystal neutron diffraction 
data is available for the former two compounds,!^^'^^^ 
information about the magnetic structure of NaCu2 02 
has thus far only been gleaned from more limited neu- 
tron powder diffractioUjl^SI nuclear magnetic resonance 
(NMR) ,[211221 and resonant x-ray diffractiorPI data. Moti- 
vated by the controversy about the origin of multiferroc- 
ity in copper-oxide chain compounds outlined above, we 
have carried out a single-crystal neutron diffraction study 
of the magnetic structure of NaCu202. We present the 
results in the framework of a rigorous symmetry classifi- 
cation of possible magnetic structures. 



II. RESULTS AND ANALYSIS 

Single crystals of NaCu2 02 in the shape of small rect- 
angular plates a few millimeters long and wide and 
~ 0.1 mm thick were grown by a self-ffux technique 
described elsewhere.'22 X-ray diffraction]^ induction- 
coupled plasma atomic emission spectroscopy,^ and 
magnetometry=231were used to fully characterize the struc- 
ture, purity, and quality of the crystals. Neutron diffrac- 
tion experiments were carried out at the Institut Laue- 
Langevin using the high-resolution powder diffractometer 
D2B, the Laue single-crystal diffractometer Vivaldi, and 
the four-circle diffractometer DIO. For the D2B measure- 
ments, we used a powder sample described elesewhere.'^S 
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FIG. 1: Powder diffraction data on NaCu202 taken on D2B 
with a neutron wavelength A = 1.596 A at T = 1.5 K. The line 
through the data is the result of a standard Rietveld refine- 
ment. The lower curve indicates the corresponding residuals. 
Tick marks indicate the positions of nuclear Bragg reflections. 



Figure 1 shows powder diffraction data taken on D2B, 
which were used to accurately refine the nuclear struc- 
ture at temperatures 1.5 and 300 K. Rietveld refinement 
in the orthorhombic space group Pnma (No. 62) with 
parameters shown in Table I yields an excellent descrip- 
tion of the data (Fig. 1). 

TABLE I: Lattice parameters and atomic positions in frac- 
tional coordinates, as derived from Rietveld refinement of 
powder diffraction data at temperatures T = 1.5 and 300 K. 
The space group is Pnma. All atoms are in Wyckoff position 
4c. 



r (K) 




1.5 


300 


a (A) 




6.2001(1) 


6.2148(1) 


fe(A) 




2.9310(1) 


2.9361(1) 


c(A) 




13.0337(2) 


13.0731(3) 


Cu(I) 


X 


0.3810(7) 


0.3811(1) 




z 


0.2549(2) 


0.2556(2) 


Cu(II) 


X 


0.8694(7) 


0.8691(8) 




z 


0.6059(2) 


0.6050(2) 


Na 


X 


0.358(1) 


0.355(1) 




z 


0.5806(4) 


0.5800(5) 


0(1) 


X 


0.4189(6) 


0.4195(7) 




z 


0.1155(3) 


0.1160(3) 


0(11) 


X 


0.3489(8) 


0.3491(1) 




z 


0.3979(3) 


0.3980(3) 


Rf 




0.054 


0.055 



In order to determine the magnetic structure of 
NaCu2 02, a set of 40 magnetic reflections was acquired 
on a single crystal at T = 1.5 K. The data were taken in 
a three-axis configuration on DIO with a neutron wave- 
length of 2.36 A. The magnetic refiections were found to 
be incommensurate with the crystal lattice and could be 
indexed as (/i, k, I) m — {h, fc, Z)jv ± k, where N stands for 
nuclear reflections. The propagation vector k ~ (0.5, ^, 0) 
with ^ = 0.228 ± 0.002 resulting from refinement of the 
data is consistent with prior neutron powder diffractiorPSl 
and resonant x-ray diffraction^'^ measurements. The tem- 
perature dependence of the magnetic Bragg intensity de- 
termined on Vivaldi (Fig. 2) showed a Neel te mpera ture 
Tjv = 12 K, again in agreement with prior work .1221231 xhe 
integrated intensities of the magnetic reflections were de- 
termined by rocking scans (inset in Fig. 2). For the 
refinement of the magnetic structure we used the mag- 
netic structure factors of 20 incommensurate magnetic 
reflections, 10 of which are unique. The overall scale fac- 
tor as well as the extinction parameter were fixed to the 
values refined from a set of nuclear reflections collected 
under the same experimental conditions. The atomic po- 
sitions and the lattice parameters were fixed to the values 
refined from the powder data (Table I). 

In order to examine the constraints set on the magnetic 
structure by the crystal symmetry, we used the represen- 
tation analysis methocP^ as implemented in the program 
Baslreps.^® The lattice symmetry (space group Pnma) 
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FIG. 2: Temperature dependence of the intensity of the mag- 
netic Bragg reflection (0.5, —0.228, —1) in NaCu202, taken on 
the Lane diffractometer Vivaldi. The inset shows a rocking 
curve taken at T = 2 K on the four-circle diffractometer DIO. 



and atomic positions derived from the structural analy- 
sis described above, as well as the propagation vector k 
determined from the positions of the magnetic reflections 
serve as input for this analysis. Further, x-ray absorption 
measurements indicate valence states of 1-1- (with a full- 
shell configuration), and 2+ (with one unpaired hole in 
the d-shell) for the Cu(I) and Cu(II) ions, respectively.!^ 
Only the latter ions therefore carry a magnetic moment. 
They are in the following positions: 



Cul 
Cu2 
Cu3 
Cu4 



(x, 1/4, z) 

{-X, 3/4, -z) + (l, 0, 1) 

{1/2 + x, 1/4, l/2-z) + (-l, 0, 1) 

(1/2 -x, 3/4, l/2 + z) + (l, 0, -1) 



with refined parameters x and z listed in Table L 

The representation analysis yields a single, two- 
dimensional irreducible representation, T 



mag : 



for the 

magnetic structure. The possible spin configurations of 
the four Cu^"*" ions in the primitive unit cell are described 
by the basis vectors (with n = 1...12) of Tmag, which 
are displayed in Table II. We denote ferro- and antifer- 
romagnetic alignment of spins along the crystallographic 
axis i by the modes /j and , respectively. The modes fx , 
tty, and fz resulting from the representation analysis im- 
ply the sign sequences (-I--I--I— (-1 1 — ), and (-I-H ) 

for the X-, y-, and z-axis components, respectively, of the 
magnetization of atoms (Cul Cu2 Cu3 Cu4). Note that 
the two-dimensionality of Tmag (which derives from the 
spin reversal along a implied by the propagation vector) 
permits a physically equivalent description in terms of 
the modes a^, /j,, and of spins in an adjacent subccU 
of the magnetic lattice (Table II). 

The observed magnetic structure factors were then 
compared to those generated by the complete set of four 
magnetic models compatible with the lattice symmetry, 
which can be specified by the possible permutations of 



TABLE II: Basis functions 'I'n of the irreducible representa- 
tion Vmag for propagation vector k = (0.5, 0.228, 0) in the 
orthorhombic space group Pnma. Only the real parts of the 
basis vectors are given, with a = cos 27r(A;j;/2) = 0.754. The 
modes ft and a; characterize ferromagnetic and antiferromag- 
netic spin configurations, respectively, along crystal axis i. 
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the real and imaginary parts, Ri and 7^, of the Fourier co- 
efficients of the magnetic moment along the crystal axes i. 
They include a sinusoidal modulation with only real com- 
ponents RxRyRz (model 1), and three elliptical helices 
with components I^RyRz, RxRylz, and RxIyRz (mod- 
els 2-4). Note that replacing R ^ I yields physically 
equivalent descriptions. In all cases, the representation 
analysis fully constrains the spin sequences of the atoms 
inside the primitive unit cell. The model proposed by 
Kobayashi and coworkers for LiCu202,'^ in which the 
phase relation of these moments was chosen arbitrarily, 
is not compatible with the lattice symmetry and was not 
considered here. The description of more complex struc- 
tures such as conical helices would require more than one 
magnetic representation. Such models are therefore also 
incompatible with our analysis, which yields a single rep- 
resentation. 

The program FullproP^ was used to perform least- 
squares refinements of the three Fourier components 
of the magnetization in the four symmetry-compatible 
structures. The outcome is shown in Table III, along with 
the average moment and the moment amplitude derived 
from these quantities and the quality-of-fit parameters 
Rf and x^. Low residuals Rp were obtained for models 
1-3. For the sine- wave modulated structure (model 1 in 
Table III) the components along the h- and c-axes were 
found to be very similar, while the component along a 
is slightly reduced. For the helical structures (models 2 
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TABLE III: Comparison of different models for the magnetic structure of NaCu2 02. The imaginary and real components of the 
Fourier coefficients of the Cu^^ magnetic moments, (7x, ly, Iz) and {Rx, Ry, Rz), as well as the resulting real part R and total 
moment M are given in units of /is. For the sinusoidal structure, the magnetic moment amplitude is listed (*). The agreement 

factors of the least-square fits arc defined as Rf = J2(\Fobs\ - \Fcaic\)/J2 \Fobs\ and = J2i\^obs — Icalc\/<y^)/(n —p), where 
n and p are the numbers of observations and refined parameters, respectively. 
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and 3 in Tabic III), wc obtained nearly the same val- 
ues of the three magnetization components as calculated 
for the sinusoidal structure. In model 2, where the real 
components are in the &c-plane and the imaginary com- 
ponent points along the a-axis, the copper moments are 
rotating in two different planes, which subtend angles of 
±45(4)° with the 6-axis and encompass the moments of 
Cul and Cu4, and Cu2 and Cu3, respectively. In model 
3, where the real components are in the a6-plane and the 
imaginary component points along the c-axis, we also ob- 
tained two different rotation planes for the moments of 
Cul and Cu3, and Cu2 and Cu4, respectively, which were 
found to subtend angles of ±50(4)° with the a-axis. In 
Table IV, the magnetic structure factors calculated for 
the latter model (which is the most likely ground state 
of NaCu202, see Section III) arc compared to the exper- 
imentally observed ones. Model 4 yields a substantially 
larger Rp and can hence be excluded. 

Figure 3 provides pictorial representations of models 2 
and 3. Both structures are elliptical helices by symmetry, 
but since the real and imaginary Fourier coefJicients are 
identical within the standard deviation of the refinement, 
the ellipticity is small. For comparison, wc also inves- 
tigated simple circular helices with polarization planes 
coincident with high-symmetry planes of the crystal lat- 
tice (models 5-7 in Table III), but these models yielded 
unsatisfactory residuals. 



TABLE IV: Comparison of the observed magnetic structure 

factors, -Fobs, of magnetic satellite reflections {h, k, 1)m and 
the structure factor F^^i^ calculated magnetic structure fac- 
tors of the spiral model 3 (see Table III). 
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FIG. 3: (Color online) Proposed magnetic structures for 
NaCu202. Model 3 yields competitive residuals in the re- 
finement of the neutron diffraction data (Table III) and is 
consistent with other experimental data on NaCu2 02 (see the 
text). The arrows indicate the sense of rotation of helices on 
different copper-oxide chains. 



III. DISCUSSION AND CONCLUSIONS 

We first discuss possible limitations of the representa- 
tion analysis that underlies the choice of magnetic struc- 
tures we have selected for comparison with the data. The 
analysis is based on the assumption that the spin Hamil- 
tonian includes terms up to bilinear order in the spin 
operators. While this is sufficient in the vast majority 
of magnetic insulators, we note that strong charge and/or 
orbital fluctuations may l ead to higher-order terms that 
require a modified analysis view of the large Mott- 
Hubbard gap of copper-oxide chain compounds^^^ and 
the large crystal-field splitting of the Cu c?-orbitals, such 
fluctuations are expected to be negligible in NaCu202. 
Since recent research has focused attention on ring- 
exchange terms in the sp in Hamiltonian of cuprate spin- 
ladder compounds J^SMl one might also consider three- 
spin interactions of moments on directly adjacent copper 
oxide chains in NaCu2 02. While no quantitative infor- 
mation on the magnitude of such interactions is available, 
they are expected to be substantially smaller than the bi- 
linear interactions and therefore have no major influence 
on the magnetic structure. An analysis of the current set 
of neutron diffraction data based on bilinear exchange in- 
teractions therefore appears to be adequate. In any case, 
an analysis based on a choice of models without recourse 
to symmetry consider at ions seems inadequate. 

Three of the four magnetic structures revealed by the 
representation analysis yield good agreement with the 
neutron data, and none of them can be singled out based 
on the diffraction study alone. We therefore discuss these 
structures in the light of data collected by other experi- 
mental probes. Amplitude-modulated states such as the 
sinusoidal structure of model 1 have been observed in 



compounds with doped edge-sharing copper oxide chains, 
which support low-energy charge fluctuationSjl^^ but as 
mentioned above, such fluctuations are strongly sup- 
pressed in Mott insulators such as NaCu202. Collinear 
incommensurate structures with weaker amplitude mod- 
ulations, such as the soliton lattice observed in spin- 
Peierls systems in high fields,'^ differ from model 1 by the 
content of higher harmonics, which could not be deter- 
mined in the present study because of the weakness of the 
corresponding higher-order Bragg reflections. However, 
the strongly anisotropic behavior of the Na NMR line- 
shape for H applied along the different crystallographic 
direction^^U jg difficult to reconcile with a collinear struc- 
ture and indicates helical order. 

The temperature dependence of the uniform magnetic 
susceptibility, Xj yields further insight into the magnetic 
order. Recent measurements on crystals from the same 
batch as the ones investigated here^^ indicate a suppres- 
sion of X for H II 6 and c, but not a, upon cooling below 
T/y. Helical states such as model 2, in which the a-axis is 
in the plane of polarization, are inconsistent with these 
data. In model 3, on the other hand, the a-axis subtends 
a larger angle with the polarization plane than both 6- 
and c-axes, in qualitative agreement with the low-field 
susceptibility data. The low symmetry of the polariza- 
tion plane in this model also explains the apparen t ab- 
sence of spin-flop transitions for fields up to 7 T.'^SMl 
Rather than a sharp spin-fiop, an external magnetic field 
along a is expected to induce a gradual rotation of the 
polarization plane towards the 5c-plane, which explains 
the good agreement of high-field NMR data with a model 
based on 6c-polarizcd spirals.^ 

The elliptical helix structure with alternating polar- 
ization planes shown in Fig. 3b is more complex than 
the simple circular, 6c-polarized helix previously identi- 
fied based on less complete powder neutron diffractiorP^ 
and NMRp2ldata. We stress, however, that these features 
are mandated by the lattice symmetry and the propaga- 
tion vector, which are accurately known for NaCu202. 
We also note the alternating sense of rotation of helices 
propagating along different copper-oxide chains (arrows 
in Fig. 3). According to Eq. 1, this implies an anti- 
ferroelectric state in which every chain generates a ferro- 
electric moment, but the macroscopic electric polariza- 
tion vanishes. This explains the absence of ferroelectric- 
ity in NaCu2 02.'^ A small reduction of the dielectric con- 
stant below T/v may be indicative of an anti-ferroelectric 
stated 

Our results also cast light on the origin of the ferro- 
electric polarization in LiCu2 02, which exhibits the same 
lattice symmetry and an incommensurate helix propa- 
gation vector of the same form, k = (0.5,^,0), as in 
NaCu202.'^ Our representation analysis therefore also 
applies to LiCu202. Since according to Eq. 1 none of the 
four magnetic states revealed by this analysis supports a 
macroscopic ferroelectric p olari zation, we conclude that 
the ferroelectricity observecF^in LiCu202 cannot be of 
intrinsic origin, and that defects generated by Li-Cu in- 
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tersubstitution must play a central role. In this respect, 
our conclusion agrees with those of Ref. [HI but disagrees 
with those of Ref. [TH We emphasize, however, that our 
findings do not invalidate models according to which heli- 
cal magnetism generates ferroelectricity.lI"2l In the frame- 
work of this scenario, substitutional defects may locally 
lift the compensation of ferroelectric moments in differ- 
ent copper-oxide chains, and a description of the resulting 
magneto-electric defect pattern may be the key to an ex- 
planation of the puzzling magnetic field dep endence of 
the ferroelectric polarization in LiCu202.'^'^ An assess- 
ment of the possible influence of defects on the ferroelec- 



tric properties of LiCuV04j^'i^ as well as other multifer- 
roics with helicoidal order,!^ is an interesting subject of 
further investigation. 



ACKNOWLEDGMENTS 

We acknowledge useful discussions with Ph. Leininger, 
G. Mclntyre, and J. Rodriguez Carvajal. We also ac- 
knowledge financial support by the DFG under grant 
SFB/TRR 80. 



^ H. Katsura, N. Nagaosa, and A.V. Balatsky, Phys. Rev. 

Lett. 95, 057205 (2005). 
^ I. A. Sergienko and E. Dagotto, Phys. Rev. B 73, 094434 

(2006) . 

^ M. Mostovoy, Phys. Rev. Lett. 96 067601 (2006). 

■* For a recent review, see K.F. Wang, J.M. Liu, and Z.F. 

Ren, Adv. Phys. 58, 321 (2009). 
^ S. Park, Y.J. Choi, C.L. Zhang, and S.-W. Cheong, Phys. 

Rev. Lett. 98, 057601 (2007). 
" Y. Naito, K. Sato, Y. Yasui, Yu. Kobayashi, Yo. 

Kobayashi, and M. Sato, J. Pliys. Soc. Jpn. 76, 023708 

(2007) . 

R.V. Pisarev, A.S. Moskvin, A. M. Kalashnikova, A.A. 

Bush, and Th. Rasing, Phys. Rev. B 74, 132509 (2006). 
® J. Malek, S.-L. Drechsler, U. Nitzsche, H. Rosner, and H. 

Eschrig, Phys. Rev. B 78, 060508(R) (2008). 
^ Y. Matiks, P. Horsch, R. K. Kremer, B. Keimer, and A.V. 

Boris, Phys. Rev. Lett. 103, 187401 (2009). 
° T. Masuda, A. Zheludev, B. Roessh, A. Bush, M. Markina, 

and A. Vasiliev, Phys. Rev. B 72, 014405 (2005). 
^ S. Seki, Y. Yamasaki, M. Soda, M. Matsuura, K. Hirota, 

and Y. Tokura, Phys. Rev. Lett. 100, 127201 (2008). 
^ Y. Yasui, K. Sato, Y. Kobayashi, and M. Sato, J. Phys. 

Soc. Jpn. 78, 084720 (2009). 
^ C. Fang, T. Datta, and J. Hu, Phys. Rev. B 79, 014107 

(2009). 

" Y. Kobayashi, K. Sato, Y. Yasui, T. Moyoshi, M. Sato, 

and K. Kakurai, J. Phys. Soc. Jpn. 78, 084721 (2009). 
^ B.J. Gibson, R. K. Kremer, A.V. Prokofiev, W. Assmus, 

and G.J. Mclntyre, Physica B 350, E253 (2004). 

Y. Yasui, Y. Naito, K. Sato, T. Moyoshi, M. Sato, and K. 

Kakurai, J. Phys. Soc. Jpn. 77, 023712 (2008). 
^ A.S. Moskvin and S.-L. Drechsler, Europhys. Lett. 81, 

57004 (2008). 

^ A.S. Moskvin and S.-L. Drechsler, Phys. Rev. B 78, 024102 

(2008) ; A.S. Moskvin, Y.D. Panov, and S.-L. Drechsler, 
ihid. 79, 104112 (2009). 

^ V. Prokofiev, I.G. Vasilyeva, V.N. Ikorskii, V.V. Malakhov, 



LP. Asanov, and W. Assmus, J. Solid State Chem. 177, 
3131 (2004). 

^° L. Capogna, M. Mayr, P. Horsch, M. Raichle, R. K. Kre- 
mer, M. Sofin, A. Maljuk, M. Jansen, and B. Keimer, Phys. 
Rev. B 71, 140402(R) (2005). 

M. Horvatic, C. Berthier, F. Tedoldi, A. Comment, M. 
Sofin, M. Jansen, and R. Stern, Prog. Theor. Phys. Suppl. 
159, 106 (2005). 

A.A. Gippius, A.S. Moskvin, and S.-L. Drechsler, Phys. 
Rev. B 77, 180403(R) (2008). 

Ph. Leininger, M. Rahlenbeck, M. Raichle, B. Bohnenbuck, 
A. Maljuk, C. T. Lin, B. Keimer, E. Weschke, E. Schierle, 
S. Seki, Y. Tokura, and J. W. Freeland, Phys. Rev. B. 81, 
085111 (2010). 

A. Maljuk, A. B. Kulakov, M. Sofin L. Capogna, J. 
Strempfer, C. T. Lina, M. Jansen, and B. Keimer, J. Cryst. 
Growth 263, 338 (2004). 

E. F. Bertaut, Acta Cryst. A 24, 217 (1968); J. Magn. 
Magn. Mater. 24, 267 (1981). 
'^^ J. Rodriguez Carvajal, Physica B 192, 55 (1993). 

B. Bohnenbuck, I. Zegkinoglou, J. Strempfer, C.S. Nelson, 
H.H. Wu, C. Schiifiler-Langeheine, M. Reehuis, E. Schierle, 
Ph. Leininger, T. Herrmannsdorfer, J.C. Lang, G. Srajer, 

C. T. Lin, and B. Keimer, Phys. Rev. Lett. 102, 037205 
(2009). 

M. Reehuis, C. Ulrich, P. Pattison, B. Ouladdiaf, M. C. 
Rheinstadter, M. Ohl, L. P. Regnauh, M. Miyasaka, Y. 
Tokura, and B. Keimer, Phys. Rev. B 73, 094440 (2006). 
S. Brehmer, H.-J. Mikeska, M. Miiller, N. Nagaosa, and S. 
Uchida, Phys. Rev. B 60, 329 (1999). 
^° A. Gofihng, U. Kuhlmann, C. Thomsen, A. Lofi'ert, C. 
Gross, and W. Assmus, Phys. Rev. B 67, 052403 (2003). 
M. Raichle, M. Reehuis, G. Andre, L. Capogna, M. Sofin, 
M. Jansen, and B. Keimer, Phys. Rev. Lett. 101, 047202 
(2008). 

V. Kiryukhin, B. Keimer, J. P. Hill, and A. VigUante, Phys. 
Rev. Lett. 76, 4608 (1996). 



